Theory of Normal State Pseudogap Behavior in FeSei_ x Te 
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The normal state of the recently discovered Iron Selenide (FeSe)-based superconductors shows a 
range of inexplicable features. Along with bad-metallic resistivity, characteristic pseudogap features 
and proximity to insulating states, reminiscent of the underdoped high-T c cuprates, mark these sys- 
tems as strongly correlated non-Fermi Liquid metals. Here, using the first-principles LDA+DMFT 
method, we show how strong multi-orbital correlation-induced orbital-selective Mott-like physics 
leads to an orthogonality catastrophe underpinning these inexplicable incoherent features. Excel- 
lent agreement with a range of resistivity and Seebeck data strongly support our proposal. We 
discuss pseudogap regime microscopically, along with implications for the nature of the instability 
at lower T, and propose that related systems could be of use in thermoelectric devices. 



The Iron Selenides (FeSe), with or without Tellurium 
(Te) substitution, are the latest addition to a rapidly 
growing list of Fc-based unconventional superconduc- 
tors These latter systems are structurally simpler 
than their 1111-oxyarsenide counterparts Q; they never- 
theless exhibit a host of very unusual physical responses. 
These are (by no means a complete list): (i) a linear- 
in-T dependence of the bad metallic normal state resis- 
tivity, pdc{T) ~ po + AT above T c , and a Te- and Cu- 
substitution induced metal-insulator transition (MIT) in 
FeSei-^Te^ Q and Fe^Cu^Se [|. (U) The NMR relax- 
ation rate is distinctly of the non-Korringa form, and the 
static (uniform, q = 0) spin susceptibility, x(T) c± T 1+n 
with < n < 1, in FeSe |5|. (Hi) An almost non- 
existent Drude peak, with dominant incoherent features 
in the optical conductivity [f|. (iv) Incoherent pseudo- 
gap like low-energy features in angle-integrated photoe- 
mission (PES) 0. Also, an ARPES study for the anti- 
ferromagnetic compositions shows that the AF ordering 
wave- vector, Qaf, is very different from that predicted 
by local-density-approximation (LDA) 

In particular, susceptibility and PES data show up the 
normal state pseudogap (PG) in FeSei-^Tea;. Several of 
the above features are reminiscent of the high-T c cuprates 
in the under- to optimally doped regime, and, taken to- 
gether, defy an interpretation in terms of the Landau 
Fermi Liquid (LFL). These PG features must manifest 
themselves in other probes as well. In fact, recent nor- 
mal state resistivity and Seebeck data also exhibit PG 
features which closely correlate with each other [the PG 
in Pdc(T), as seen in (dp/dT) and (d 2 p/dT 2 ) near the PG 
scale T* = Tpq, correlates well with the related anomaly 
in S(T) at the same T = T*] Further, at low T, 

the thermopower does not show the S(T) ~ aT behav- 
ior; the T dependence is distinctly slower than linear just 
above T c . Along with Pdc(T) ~ AT, this confirms that 
the normal state in FeSei-^Te^ is not a LFL metal. In 
addition, S(T) at high-T flattens out 0, being well de- 
scribed by a Heikes-like law; this is exactly the behavior 
expected from a Hubbard-like model at high T. Finally, 



S(T, x) data 0] reveal an approximate isosbectic point at 
T ~ 100 K. This is strong evidence for a crossover scale, 
associated with effectively localized (Heikes law in S(T)) 
to incoherent metallic conduction, in FeSei-^Te^. 

Taken together, the above also imply that the 
FeSei-^Te^ are strongly correlated materials proximate 
to a Mott insulator. The microscopic electronic cor- 
relation processes underlying the emergence of an in- 
coherent pseudogapped metallic phase above T c (Tjv) 
in FeSei-ajTe^, however, await a consistent theoretical 
understanding. Like the 1111-systems, Fe-selenides are 
multi-orbital (MO) systems with all five d bands cross- 
ing the Fermi energy (Ep) in band (LDA) calculations. 
Given the global non-LFL features found experimen- 
tally, however, non-trivial extensions of the LDA to ac- 
count for dynamical correlation-driven incoherence and 
breakdown of LFL theory are mandatory. Here, using 
the state-of-the-art LDA+dynamical-mean-field-theory 
(LDA+DMFT), already used with good success [ll| for a 
range of correlated systems, we study this pseudogapped 
incoherent metallic state in FeSei-^Tcj; in detail. 

In earlier work, we have shown that the incoher- 
ent metal [l2} as well as the Te-doping induced metal- 
insulator transition (MIT) (l3| can be quantitatively un- 
derstood using LDA+DMFT with sizable d band MO 
correlations. Here, we extend these to characterize the 
PG phase in depth. Specifically, we show how the See- 
beck co-efficient (thermopower), S(T) is also quantita- 
tively described, and correlate the specific PG anomalies 
in both S(T) and Pdc(T) with doping, x. Very good 
quantitative accord with extant data reinforces the ba- 
sic hypothesis about crucial role of strong correlations. 
To our best knowledge, this is a first attempt to uncover 
such features in an incoherent metal within LDA+DMFT 
for a real system. Finally, we show how the itinerant- 
localized duality (Mottness) makes for a high thermo- 
electric figure-of-merit (ZT) in such systems at low T, 
and suggest possibilities for thermoelectric applications 
using related systems. 

Since the basic LDA+DMFT formulation has been 
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used for a variety of real systems with good quantita- 
tive success [ll|, we do not present it here, but directly 
describe the electronic transport within LDA+DMFT 
scheme. The Seebeck coefficient, like the conductivity, is 
exactly computable within DMFT using the fully renor- 
malized LDA+DMFT propagators. We generalized the 
DMFT result [HI, Ell to tne hve-band case relevant for 
FeSei-^Tea;. This involves the following steps: 

We begin with the general expressions for Seebeck 
and thermal conductivity, which, respectively, measures 
the mixed electrical-thermal correlations [<S(T)] and the 
heat-current correlations [k(T)} at finite T: 
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with A n (T) = ^ /+~ dw^(w)[-/'(w)](w - A*)" and 
r x (u) = -^^Tr[^, a (k)p a (k, W ) U:c , fc (k)p h (k, W )] . 

k 



(3) 

is the LDA+DMFT 
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Here, p„(k, = ——Im — ^ / , 

spectral function, a = (xy, xz, yz, x 2 — y 2 , 3z 2 — r 2 ) rep- 
resent the five d orbitals, p is the chemical potential, 
u x ,a(k) is the group velocity, /(w) is the Fermi-Dirac 
function and V the sample volume. In DMFT, the A n (T) 
are convertible to integrals over the unperturbed DOS. 
The only approximation made here is to ignore the k 
dependence of v x a (k), i.e, v x a (k) — > v x>a = v. In an 
incoherent metal, such as we have here, this is justified, 
since, between successive hops, a carrier in an incoherent 
state does not exist long enough in a given k eigenstate. 
Another related reason is the bad-metallic resistivity, im- 
plying UfI — 0(1) and leading to the same conclusion. 
With this simplification, Eq. ([3|) for the intra-band con- 
tribution reads (j) xx (w) = jrJ2akPl( k ^)- 

We now describe our results. In Fig. [TJ we show the 
resistivity, PdciT) = 1/Aq(T) in FeSe^^Te^, derived ear- 
lier [l3j]. It clearly shows the "S-like" shape character- 
istic of a PG metal. At low T, p dc (T) ~ AT and the 
LFL-like T 2 form is never observed. With increasing x, 
Pdc becomes more bad-metallic, smoothly going over to 
an insulator-like form up to very low T ~ 2.0 K. In re- 
ality, at much higher T, either antiferromagnetic (AF) 
or unconventional superconducting (USC) transition will 
cut-off this extremely low crossover to a metal. For non 
AF/USC ordered systems, even a minute amount of non- 
magnetic disorder on FeSe layers (e.g, Cu 0] will imme- 
diately destroy metallicity, leading to a metal-insulator 
transition, as indeed seen in Fei-zCuzSe for x « 1. Fo- 
cussing on x — to make contact with experiment [loj . 
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FIG. 1: (Color online) Resistivity versus temperature (lower 
panel) and its derivatives (upper panel) showing their evolu- 
tion with increasing x (n t — 6 + x) within LDA+DMFT for 
fixed U — 4.0 eV and Jh =0.7 eV. The pseudogap scale is 
Tpg 150 K (see text). 



we show (dpdc/dT) and (d 2 pd c /dT 2 ) in the top pan- 
els of Fig. [TJ The PG scale, T* is the temperature at 
which (dpdc/dT) shows a maximum, while (d 2 pd c /dT 2 ) 
crosses zero [nj. The Seebeck co-efficient, S(T), should 
also show an anomalous feature, i.e, a minimum, around 

t* HIES- 

Very good quantitative agreement with extant data [^, 
lioj is also clearly visible in S(T, x), Fig.O Our results are 
obviously valid only for T > T c , and, at very low T, cor- 
respond to what would be observed in a non-SC system 
in absence of disorder. A number of characteristic fea- 
tures, intimately correlated with those in the resistivity 
derivatives, are noteworthy: (i) S(T) is approximately T- 
independent at T > T m j n (x), in nice accord with Hcikcs 
law. Here, it arises due to thermal transport involving 
Mott localized carriers in our five-orbital Hubbard model. 
(m) In the relevant range of < x < 0.6, we find that 
S(T) smoothly decreases, going through a minimum at 
another scale T*(x). Interestingly, the deviation from 
Heikes law sets in around Tpg — 150 K. Further, this 
scale coincides with the PG scale extracted from resis- 
tivity derivatives (see Fig. [JJ, pinning down the temper- 
ature below which a low energy PG opens up in trans- 
port. (Hi) S(T,x) crosses zero twice at Ti ow (x),Thi g h(x) 
(labelled T an , T* in Ref. [HJ), again m precise accord 
with data, (iv) Even at low T ~ 10 K, S(T) does not 
recover the linear-in-T form: an examination shows a 
T v law with rj < 1, suggesting that collective bosonic 
fluctuations dominate the low-T thermopower. (v) Also 
the approximate isosbectic point in S(T, x) as x is varied 
is reproduced by LDA+DMFT. (vi) Interestingly, S(T) 
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only shows two sign changes for metallic (above T c ) com- 
positions. For large x > 0.5, appearance of insulator-like 
behavior above T c [l3[ goes hand-in-hand with only one 
sign change at high T* ~ 150 K in S(T), also consistent 
with findings |9|. 
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FIG. 2: (Color online) Seebeck coefficient (top left) and 
thermal conductivity (top right) for FeSei-zTex within 
LDA+DMFT. Notice the S(T) =const. form at high-T, de- 
creasing around Tpg and showing a second crossin g at much 
lower Tan, all in excellent agreement with data [9l llOfl. The 
thermoelectric figure of merit, ZT, is shown in the lower 
panel. Notice the high value of ZT for x — 0.6. 



fall-off of the DOS at energies above the (small) PG fea- 
ture around Ep instead of a LFL peak at Ep is seen for 
the d XZ) y Z)X 2_ y 2 DOS. Finally, the low-energ kink feature 
clearly seen for x — (solid curve) is smoothened with x. 
Based on the agreement with transport above, we pre- 
dict that PES (ARPES) measurements on FeSei_ x Te x 
as a function of x should show all these features in the 
spectral function with x. 
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FIG. 3: (Color online) LDA+DMFT spectral functions for 
FeSei-^Te^. Clear pseudogap stabilization with increasing x 
going hand-in-hand with large scale SWT is apparent. Power- 
law fall-offs in the d xzyzx 2_ y 2 spectra above the low-energy 
PG, associated with orbital-selective incoherence, are also vis- 
ible. 



Remarkably, the full set of our theoretical Pdc(T, x) 
and S(T, x) results are in very good agreement with a 
range of extant results obtained by two groups [j| Hoj . 
Our work is the first detailed theoretical study of elec- 
trical and thermal transport in Fe-based superconduc- 
tors. These remarkable results clearly call for a deeper 
microscopic rationalization. Since vertex corrections to 
conductivities rigorously vanish in DMFT [l6[ , these fea- 
tures must be intimately linked to the detailed evolution 
of the LDA+DMFT spectral functions with T and x. In 
Fig. H we show the LDA+DMFT orbital-reolved DOS, 
p a (u), for the doping values above. Clear signatures of 
normal state incoherent metal behavior, along with a 
low-energy PG, are visible in the DOS, providing a clear 
link between PG features in transport and spectral re- 
sponses. In particular, we observe important changes in 
the LDA+DMFT DOS with x, and these must be corre- 
lated with features in transport. The d XZtVZ DOS always 
has a deep PG, while the d xy DOS rapidly loses the small 
LFL coherence with x. The d x 2_ y 2 DOS shows a progres- 
sively deeper PG feature with x, while the reverse occurs 
for the ds z 2_ r 2 DOS. These go hand-in-hand with large 
changes in spectral weight transfer (SWT), a character- 
istic feature of Mottness Further, clear power- law 



Also the anomalies in S(T) and (d n pd c /dT n ) with 
n = 1,2 at the much lower T an 1 1711 . observed for the 
first time in a correlated system [10(, are recovered too 
within our LDA+DMFT. In a band-LFL picture, this 
behavior would be associated with hole conduction, but 
such a connection is tenuous in the incoherent metal we 
find. In real FeSei-^Te^, USC and/or AF actually cuts 
off this behavior, leading to S(T — > 0) = 0. Looking 
closer at the T — > limit within our uniform, symmetry- 
unbroken, state calculation, we find that S(T — )• 0) ~ T n 
with 77 < 1. Correspondingly, pd c {T) ~ AT up to lowest 
T. Neither p dc {T) ~ T nor S{T) ~ T v with 77 < 1 at 
low T are interpretable in terms of transport associated 
with LFL quasiparticles; rather, they suggest transport 
in terms of incoherent, collective excitations in the non- 
LFL metal, //the above features are related to the mech- 
anism of superconductivity (SC) [l(J, our analysis has 
far-reaching implications for the SC-instability. Specifi- 
cally, our results now imply that SC must arise directly as 
an instability of a normal non-LFL metal without long- 
lived fermionic quasiparticles, and that the d xz yz x 2_ y 2 
bands will play a central role in pairing, with the rest 
playing a secondary role, as in a proximity effect sce- 
nario. Recall that, in our earlier LDA+DMFT study 
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of FeSe 12|, [13J , the non-LFL metal arises from orbital 



selective incoherence having its origin in the interplay 
between strongly orbital-dependent hopping (as in LDA) 
and strong dynamical inter-orbital correlations (DMFT); 
the non-LFL features are assigned to the generation of an 
Anderson orthogonality catastrophe in the local problem 
of DMFT [jjj]- Once this occurs, the spectral functions 
will exhibit a low energy PG with power-law like tails 
at higher energy, precisely as we find in the DMFT. In 
such a metal, the transport properties are dominated by 
collective, multi-particle bosonic excitations flij ]. rather 
than by LFL quasiparticles. Our LDA+DMFT results 
are consistent with this incoherence scenario based on 
Mottness, and excellent agreement with data from a num- 
ber of groups, as shown above, strongly supports our con- 
tention. 

Therefore, we extend our analysis to propose that re- 
lated materials have the potentiality to be of use as good 
thermoelectrics, at least at low T. To this end, we also 
need the thermal conductivity, k(T), also computable ex- 
actly in DMFT in terms of the LDA+DMFT propagators. 
In Fig. [21 we show the electronic part of the thermal con- 
ductivity [k(T)}, along with the thermoelectric figure of 
merit, ZT = 2^2. We find jfe(T) ~ T 2 at low T, indi- 
cating that bosonic contributions dominate the electronic 
contribution to k(T). We suggest that measuring k(T) 
at low T will bare this contribution. At higher T > 10 K, 
additional phononic contributions will increase k(T) and 
decrease the electrical conductivity cr(T) = Aq(T), reduc- 
ing ZT. Thus, given the weak electron-phonon coupling 
in the Fe-based SC [20j . we expect a small-to- moderate 
reduction in ZT at low T. Nevertheless, the high value of 
ZT, albeit at low T, is encouraging, and ways of improv- 
ing ZT at higher T, based on appropriate engineering 
are an attractive possibility. In particular, consideration 
of more geometrically frustrated members of the family 
of Fe arsenides and selenides could yield an appreciable 
S{T) ~ T at high T arising from large spin and orbital 
fluctuations, as in NaCo204 [2lj, and increasing ZT. 

Given an incoherent non-LFL metal with low energy 
PG as above, the low T instabilities cannot involve usual 
particle-hole (p-h) (AF) or particle-particle (p-p) (USC) 
pairings in a LFL, since there are no coherent LFL quasi- 
particles in the non-LFL state in the first place. We be- 
lieve that these instabilities will then arise by studying 
the dominant two-particle instabilities of such a metal. 
As in coupled D = 1 Luttinger liquids [jjj . in a regime 
where one-particle inter-orbital mixing is irrelevant, gen- 
eration of effective two-particle (pair hopping) terms in 
both p-h and p-p sectors at two-loop level in a renormal- 
ization group treatment around the impurity (DMFT) 



limit lead to ordered states emerging directly from the 
incoherent metal 23[ . We leave the details of the mecha- 
nism of USC from an incoherent metal derived above to 
the future. 

In conclusion, we have studied the detailed nature of 
the incoherent metal with a low energy pseudogap using 
the first-principles LDA+DMFT. The anomalous charac- 
ter of the transport properties is interpreted in terms of 
electrical and heat transport by collective, multi-particle 
excitations. These arise from an Anderson orthogonal- 
ity catastrophe having its origin in the selective-Mott 
localization in LDA+DMFT. Excellent semiquantitative 
agreement with a whole host of transport (pd c (T) and 
S(T)) data with x puts our mechanism on solid ground. 
Our study puts strong theoretical constraints on how the 
ordered states (AF/USC) arise directly from the incoher- 
ent metal, and shows the crucial importance of Mottness 
in this context. 
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